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CLAIMS 

1. A method of determining subsidence in a reservoir, said subsidence resulting from 
a continual withdrawal of hydrocarbon deposits and other fluids from said reservoir, 
comprising the steps of: 

(a) estimating rock displacement parameters u, v, and w representing rock movement 
in the x, y, and z directions; 

(b) determining £^ (*x,y,z elongation strains*) and J^y^y.^^ ('shear strains') from the 
rock displacement parameters u, v, and w; 

(c) determining o^, ('elastic normal rock stress in x,y,z directions*) and ij^^y^^^^ 
('elastic shear stress') from the e^ ^^ and 

(d) solving a set of rock momentum balance differential equations from the o^ y , and 

the ^jty^^ and determining if any residuals exist, the rock displacement parameters u, 

V, and w representing accurate rock displacement parameters for said reservoir when 
the residuals are substantially equal to zero; and 

(e) determining said subsidence in said reservoir from said accurate rock displacement 
parameters determined from step (d), by solving a final set of failure criterion 
equations, which comprise the residuals and any derivatives, to determine a set of 
rock plastic displacements, said rock plastic displacements forming a part of said rock 
displacement parameters u, v, and w. 
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2. The method of claim 1, wherein the detemiining step (b) comprises the step of: 

determining said b^^ ^ (the *x,y,z elongation strains') and said >^ (the 'shear 

strains') from said rock displacement parameters u, v, and w by using the following 
equations: 
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3. The method of claim 2, wherein the detennining step (c) comprises the step of: 

determining said a^^ ^ (the 'elastic normal rock stress in x,y,z directions') and said 
'^.ry.yz^ (the 'elastic shear stress') from said and „ by using the following 
equations: 

= 2Gf , + a(£, + + f J 
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4. The method of claim 3, wherein the solving step (d) comprises the steps of: 

solving said set of rock momentum balance differential equations from said a, ^. , and 
said „ „ , 

v/herein said rock momentum balance differential equations have a left-hand side on a 
left side of an equal sign and a right-hand side on a right side of said equal sign, said 
rock momentum balance differential equations including. 



— ^ + — ^ + — ^ + F +P =0 
3;c dy dz 



da dr dt 

+ + ^ + F +P =0 
oz dx dy 



said rock displacement parameters u, v, and w estimated in step (a) representing said 
accurate rock displacement parameters for said reservoir when a left-hand side of said 
rock momentum balance differential equations substantially equal a right hand side of 
said rock momentum balance differential equations, 

5. The method of claim 4, wherem the solving step (e) comprises the steps of using 
the said and said "i^^^y^^^ to determine whether a failure criterion has been 

exceeded in one of a set of staggered grids. 
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6. The method of claim 5, wherein, when said failure criterion has been exceeded in 
said one of said set of staggered grids, the solving step (e) further comprises the step 
of; 

setting up a final set of equations, each equation being a failure criterion 
approximately equal to zero, and 

determining a set of parameters up, vp, wp which are plastic displacements resulting 
from a failure of the rock, said plastic displacements forming a part of said subsidence 
and being determined from said'accurate rock displacement parameters for said 
reservoir. 

7. A program storage device readable by a machine, tangibly embodying a program 
of instructions executable by the machine to perform method steps for determining 
subsidence in a reservoir, said subsidence resulting from a continual withdrawal of 
hydrocarbon deposits and other fluids from said reservoir, said method steps 
comprising: 

(a) estimating rock displacement parameters u, v, and w representing rock movement 
in the x, y, and z directions; 

(b) determining y , Cx,y,z elongation strains') and Jj^^y^^^ ('shear strains') from the 
rock displacement parameters u, v, and w; 

(c) determining o^ ^ ^ ('elastic normal rock stress in x,y,z directions') and 'ij^.^y^^^^ 
('elastic shear stress') from the a^ y ^ and 7^,3^,5^ ^ 

(d) solving a set of rock momentum balance differential equations from the o^^y^ and 
the "ixy^yi^^ and determining if any residuals exist, the rock displacement parameters u. 
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V, and w representing accurate rock displacement parameters for said reservoir when 
the residuals are substantially equal to zero; and 



(e) detennining said subsidence in said reservoir from said accurate rock displacement 
parameters detemiined from step (d), by solving a final set of failure criterion 
equations, which comprise the residuals and any derivatives, to determine a set of 
rock plastic displacements, said rock plastic displacements forming a part of said rock 
displacement parameters u, v, and w. 

8. The program storage device of claim 7, wherein the determining step (b) 
comprises the step of: 

detennining said s^ y ^ (the 'x,y,z elongation strains') and said (the 'shear 

strains') from said rock displacement parameters u, v, and w by using the following 
equations: 



dv 



dw 
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dw du 
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9. The program storage device of claim 8, wherein the determining step (c) comprises 
the step of: 

determining said o^ ^ ^ (the 'elastic normal rock stress in x,y,z directions') and said 
^xs,yi.xz ^elastic shear stress') from said and 7^,^^ by using the following 

equations: 

10. The program storage device of claim 9, wherein the solving step (d) comprises 
the steps of: 

solving said set of rock momentum balance differential equations from said ^ ^ and 
said t 

wherein said rock momentum balance differential equations have a left-hand side on a 
left side of an equal sign and a right-hand side on a right side of said equal sign, said 
rock momentum balance differential equations including. 



^CTx ^^yx _ 

By dx dz ^"^y^^y ^ 
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said rock displacement parameters u, v, and w estimated in step (a) representing said 
accurate rock displacement parameters for said reservoir when a left-hand side of said 
rock momentum balance differential equations substantially equal a right hand side of 
said rock momentum balance differential equations, 

11. The program storage device of claim 10, wherein the solving step (e) comprises 
the steps of using the said o^ ^ ^ and said ^^^y^^^z determine whether a failure 
criterion has been exceeded in one of a set of staggered grids. 

12. The program storage device of claim 11, wherein, when said failure criterion has 
been exceeded in said one of s?dd set of staggered grids, the solving step (e) further 
comprises the step of: 

setting up a final set of equations, each equation being a failure criterion 
approximately equal to zero, and 

determining a set of parameters up, vp, wp which are plastic displacements resulting 
from a failure of the rock, said plastic displacements forming a part of said subsidence 
and being determined from said accurate rock displacement parameters for said 
reservoir. 
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